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Motivated by the recent evidence for direct CP-violation in D° — > h + h~ decays, we provide 
an exhaustive study of both Cabibbo-favored and Cabibbo-suppressed (singly and doubly) D° — > 
K^h^t^ £~ decays. In particular, we study the Dalitz plot for the long-distance contributions in 



the (m H ,m hh ) parameter space. We find that near-resonant effects, i.e., D 



V{h+hz)l + t 



with V — p,K*,<f>, are sizeable and even dominant (over Bremsstrahlung) for the n + \i~ decay 
modes, bringing the branching ratios close to the LHCb reach. We also provide a detailed study 
of the angular asymmetries for such decays and identify signatures for new physics detection. In 
particular, new physics signals can be neatly isolated in asymmetries involving the semileptonic 
penguin Qio, where for typical new physics scenarios the effects can be as sizeable as 0(10%) for 
the doubly Cabibbo-suppressed modes. 



PACS numbers: 13.20.Fc, 12.39.St, 11.30.Er 

I. INTRODUCTION 

Processes involving flavor changing neutral cur- 
rents (FCNC) are loop-suppressed in the standard 
model and therefore are especially suited as probes 
of new physics. In D decays the suppression is 
even more accentuated that in B or kaon decays: 
in charm physics FCNC involve down-type quarks 
and as a result the GIM mechanism is more efficient. 
Since new physics does not have to be subject to the 
same GIM suppression, charm decays are in princi- 
ple an ideal arena to test physics beyond the stan- 
dard model. The situation is of course not so simple: 
light quarks carry the bulk contribution to the de- 
cays, which means that they are long-distance dom- 
inated. The resulting hadronic uncertainties, due to 
their nonperturbative nature, are very difficult to es- 
timate and overshadow short-distance effects, mak- 
ing their detection rather challenging. The situation 
gets worse because charm physics does not seem to 
accept an effective field theory description and one 
has to resort, for instance, to lattice computations 
or hadronic models. 

Quite recently the LHCb [1] and CDF @, U col- 
laborations reported convincing evidence for direct 
CP violation in D° —> ir + Tr~ , K + K~ decays. Specif- 
ically, they found that the CP asymmetry Aacp = 
gives a nonzero value, that averaged with 



a CP ~ 

the previous results of the B-factories d, [j| gives 

Aacp = (-0.65 ±0.18)% (1) 

A standard model interpretation of the previous re- 
sults seems hard to accomodate (see, e.g. @,Q), even 
allowing for generous uncertainties in the estimation 
of hadronic matrix elements (unless large penguin 



contributions are at work Q). 

If new physics is the explanation behind Eq. ([1} , 
then one should scrutinize other decay modes in 
search of similar large effects. For instance, large 
CP-violating effects for radiative D° — > V-f decays 
would naturally point out at sizeable electromag- 
netic penguins 110], which is a rather common fea- 
ture of many extensions of the standard model. 

In this paper we will study the semileptonic 4- 
body decays D° —> h\hl 1 l Jr l~ , (I = e, p) in all its 
variants: Cabibbo allowed (K~ir + ), singly Cabibbo- 
suppressed (7r + 7r~; K + K~) and doubly Cabibbo- 
suppressed (K + ir~). Such a study is rather timely: 
LHCb has recently reported the potential for reach- 
ing branching ratios of 10~ 6 in the dimuon de- 
cays at the 3cr level [Tlj . improving previous up- 
per bounds [l2l [l3[ by one order of magnitude. 
BESIII should also be able to study 4-body 
semileptonic D° decays in the near future. 

On the theoretical side, there are a number of mo- 
tivations to study these rare decays: 

• The angular structure of a 4-body decay allows 
to define a variety of differential distributions. 
The associated Dalitz plots become essential 
tools for detailed tests of both the standard 
model and new physics. Of special importance 
are the different angular asymmetries that one 
can construct, which allow for a clean separa- 
tion of short and long-distance effects. 

• Access to short-distance physics is not lim- 
ited to the charge asymmetry. In particular, 
some observables offer the possibility of disen- 
tangling new physics contributions in observ- 
ables with tiny standard model backgrounds, 
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like forward-backward asymmetries. This is 
especially interesting for the doubly Cabibbo- 
suppressed mode, which is particularly sensi- 
tive to new physics. 

In a nutshell, the penalty of small branching frac- 
tions one naturally pays in 4-body decays as opposed 
to 2 or 3-body decays is overly compensated by the 
diversity (and the size) of the asymmetries one can 
build. 

The only reference for D° — > hth^£ + £~ in the 
literature is the recent work in |15| . where generic 
estimates are given for the branching ratios. In 
this work we will refine these estimates for the eight 
channels under study, detailing the differential dis- 
tributions of the different long-distance contribu- 
tions (bremsstrahlung and both magnetic and elec- 



tric hadronic components) in the (mj 



l hh. 



plane. 



Our results for the branching ratios are substan- 
tially larger than previously estimated in [l5[ : D° — > 
K + K~£ + £~ hovers around 10" 7 , while we expect 
D° tt + 7t-£ + £- to be within the reach of LHCb. 
The Cabibbo-allowed mode D° ->■ K-tt+£+£- 
is predicted in the high 10 -6 while the doubly 
Cabibbo-suppressed D° — > K + ir~ £ + £~ is estimated 
at 1CT 8 . 

We will then proceed to study angular asymme- 
tries that isolate short-distance effects. We will 
concentrate on two asymmetries that, due to tiny 
standard model backgrounds, are clean tests of new 
physics. First we will consider a T-odd asymme- 
try Ac/,, resulting from the interference between the 
electric and magnetic hadronic pieces, where is 
the angle between the dihadron and dilepton planes. 
This asymmetry involves four-quark operators and 
can be parametrized in terms of a weak phase Sw- 
Next we will consider the forward-backward asym- 
metry for the dilepton pair, Aps, which involves 
the semileptonic penguin operator Q\q. While no 
generic prediction of their magnitude can be given 
(without resorting to models of new physics), we 
show that in both cases the signal is concentrated 
around the resonant region, i.e., along the line de- 



fined by mj 



(0.5 - 1) GeV 2 and m 
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where H = p,K*,<fi for hh = tt + tt, K~it + , K + K~, 
respectively. Taking some reference values for Sw 
and Cio we show that A^ ~ (1 — 8)%, where the 
D° — > ir + ir~£ + £~ modes are the most favored, 
while Afb can reach 10% for the doubly Cabibbo- 
suppressed decay D° —> K + tt~ n + p~ . 

We will organize this paper as follows: in 
Section |TT] we describe long distances, both 
Bremsstrahlung and hadronic contributions, and 
discuss their properties for the different decay modes 
in the (w^m 2 ^) Dalitz plot. In Section ITO1 we re- 
view the short-distance effects to D° —> 2h2£ within 
the standard model. Section [IV] is devoted to an- 



D n (P) 




FIG. 1: Photon-mediated D° — > K^h'^t t t~ decay with 
our kinematic conventions. The blob represents the 
hadronic tensor H„ . 



gular asymmetries, where we study A$ and Afb as 
two examples of clean tests of new physics. Conclu- 
sions are given in Section [V] while technical details 
are collected in three Appendices. 



II. LONG-DISTANCE HADRONIC 
CONTRIBUTIONS 

As mentioned in the Introduction, D° —> 
hfh2£ + £~ are largely dominated by long-distance 
effects. The bulk of the decay width thus comes 
from light down-type quarks running into one-loop 
diagrams. At low energies, when a s ~ 1 and dy- 
namics become nonperturbativc, the picture that 
applies is shown in Fig. [TJ where the blob collects 
hadronized strange and down quarks. At the same 
time, the dilepton pair creation is dominated by pho- 
ton exchange. The amplitude for D° — > hfh^^f* ~ > 
h^h^ £ + £~ can be parametrized as 

M LD = L»(k + ,k^)H^ Pl ,p 2 ,q) (2) 

where L M is the leptonic current 

L*(fc+,fc_) = 4[fi(fc-)7'M*+)] (3) 

and Hfj, is the hadronic vector, which can be written 
in terms of three form factors Ff. 

H»( Pl ,p 2 ,q) = F lP ^+F 2 p^ + F 3 e^p lu p 2a q^ (4) 



In terms of the momenta pi, 2 , q = k + + k_ and 
Q = k + — fc_ , the squared amplitude can be cast as 



spins 



3 3 

J £\F l \ 2 T u + 2ReY,(F*F J )T l3 

i i<j 

(5) 
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where (i, j — 1,2) 



implies that 



Tij = q ■ p t q -Pj-Q- PiQ ■ Pj - q 2 Pi ■ Pj] 
T i3 = (Q ■ Pi)e^x p PiP2q X Q p ; 

T33 = 4mj [{m 2 hl m 2 h2 ~ (p! ■ p 2 ) 2 )q 2 - m 2 hl (q ■ p 2 ) 2 

- TO L(<7 -Pi) 2 + 2 Pi ■ P2q ■ piq ■ P2] 

+ (Q-Pi) 2 [(q-P2) 2 -q 2 m 2 h2 ] 
+ {Q-P2) 2 [{q-Pl) 2 ~q 2 m 2 hl ] 

+ 2Q ■ piQ ■ p 2 {q 2 pi ■ P2 - q ■ piq ■ P2) (6) 

Quite generally, there are three dynamically dis- 
tinct long-distance contributions: (i) internal 
Bremsstrahlung, i.e., QED radiation of photons 
away from the weak D° — > hfh 2 vertex; (ii) direct 
emission, where strong and weak effects get com- 
bined and resonance exchange is typically assumed 
to be the main contribution; and (iii) the charge ra- 
dius contribution, where the photon gets radiated as 
a result of D° — D° mixing. 

In the following subsections we will discuss, 
in turn, the Bremsstrahlung and direct emission 
contributions, giving expressions for the angular- 
integrated decay rates and branching ratios. The 
results will be shown in Dalitz plots for d 2 T/dq 2 dp 2 . 

In this paper we will not discuss the charge radius 
contribution since it bears no effect on the differ- 
ent angular asymmetries. Its contribution is limited 
to the decay width and subject to large uncertain- 
ties [Hj]. With this in mind, our results for the decay 
width should thus be considered as a lower bound. 
However, it's hard to imagine that the charge radius 
can increase the decay width dramatically, so for 
the purposes of this paper omitting its contribution 
should be a reasonable approximation. 



A. Bremsstrahlung 

The radiation of photons away from the weak 
D° —> h\h 2 vertex is a genuine infrared effect that 
can be computed with Low's theorem (l6| : 



M b (D° -> h+h 2l ) = 2eM(D° -> h+h 2 ) 
Pi-e p 2 ■ e 



Ff > = (-1; 



2q ■ pj + q 2 



M 



(8) 



2pi ■ q ■ 



2p 2 ■ q + q 2 



(7) 



Since the photon emission is factored out, the only 
ingredient that one needs is \A4(D° — > hth 2 )\, 
which can be extracted from experiment [6] for 
the different decay modes. Henceforth we will as- 
sume that there are no CP-violating phases in the 
D° — > h\h 2 amplitudes. 

In terms of the form factors entering H^, Eq. (J7J) 



Inserting the previous expression in Eq. ([5|), one 
can easily build the differential decay rate. The kine- 
matics of 4-body decays requires five variables. In 
this paper we will use the Cabibbo-Maksymowicz 
set of variables, Xj — (q 2 ,p 2 ,9h,0£,4>), whose con- 
crete definitions are given in Appendix [XJ In terms 
of these variables, one can ex pres s the differential 
decay rate quite generically as [13, EH 



d 5 T 
dxdy 



Ai(x) + A 2 {x)s 2 + A3(x)s 2 c^ + Ai{x)s 2 tc^, 

+ A 5 (x)s e c < j, + A e (x)c e + A7(x)s£S ( p 

+ A S (x)s 2 £S ( p + Ag(x)s'js2 < i, (9) 



where the vectors x = (q 2 ,p 2 , cos Oh) and y = 
(cos Qg, <t>) split the dynamical variables (entering the 
form factors At(x)) from the pure kincmatical an- 
gular distribution described by y. In the following 
we will integrate the full angular dependence above 
and concentrate on the quantity d 2 T/dq 2 dp 2 . De- 
tails thereof, including the analytical formulae, can 
be found in Appendix [Cl 

In Fig. |3] we show the Dalitz plot in the 
(q 2 ,p 2 ) plane for the different decay modes. Since 
Bremsstrahlung is an infrared effect driven by pho- 
ton emission, one expects a significant contribu- 
tion only in the low-g 2 region, with a sharp in- 
crease close to the dilepton threshold, whose maxi- 
mum is reached in the region of large hadron recoil 
(p 2 ~ {mo — 2mg) 2 ). In Table fl] one can read off 
the resulting branching ratios for the different de- 
cay modes. For dimuon decays, branching ratios 
span between 10 -7 — 10~ 10 , where the differences 
are solely due to the hierarchy between the hadronic 
branching ratios Bi(D° —> hih 2 ). For electron- 
positron decays there is an increase of two orders 
of magnitude per channel, with branching ratios in 
the window 10~ 5 — 1CP 8 . This increase precisely il- 
lustrates the strength of Bremsstrahlung at low-q 2 : 
electron-positron decays have a lower q 2 -threshold 
and therefore probe infrared physics deeper. Notice 
that our results for the Bremsstrahlung are substan- 
tially larger than the rough estimates given in [l5j 
by 2-3 orders of magnitude. 

Finally, we want to emphasize that the previous 
results for the Bremsstrahlung contribution only rely 
on Low's theorem (i.e., QED) and experimental in- 
put for D° — > h\h 2 . Therefore, the results in the first 
column of Table U should be seen as a solid lower es- 
timate of the total branching ratios for the different 
decay channels. 
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D 


ecay mode 


Bremsstrahlung 


Direct emission (E) 


Direct emission (M) 


D - 


T 1 + + — 


9.9 • 10~ 6 


6.2 ■ 1CT 6 


4.8 • 10~ 7 


D° - 


-¥ ir + iv~e + e~ 


5.3 • 1(T 7 


1.3 • 1(T 6 


1.3 • 1(T 7 


D°- 


> K+R-e+e- 


5.4 ■ 1(T 7 


1.1 ■ 1(T 7 


5.0- ltr 9 


D°~ 


->■ K + n~e + e~ 


3.7- 10" 8 


1.7 ■ 10" 8 


1.3 • 10" 9 


D° - 


¥ K~TT + ^l + fl~ 


8.6- 1(T 8 


6.2 ■ 10" 6 


4.8 • 10" 7 


D° - 




5.6- 1CT 9 


1.3 • 1(T 6 


1.3- 1(T 7 


D° - 




3.3 • 1(T 9 


1.1 • 1(T 7 


5.0- icn 9 


D°- 


> K + -k~ n + n~ 


3.3 • 10" 10 


1.7 • 10" 8 


1.3 • 10" 9 



TABLE I: Long-distance contributions to the branching ratio for the different decay modes. 



B. Resonant contributions 

Besides Bremsstrahlung, there are also long- 
distance contributions associated with hadronic ef- 
fects, the most important of which are the near- 
resonant regions in both the dihadron and dilepton 
sectors. Depending on the strength of the resonant 
effect, these contributions have the potential of over- 
coming the Bremsstrahlung contribution discussed 
in the previous Section. Similar contributions have 
been studied in the B -> K*(Kn)l + l~ decay pj|. 
The situation is in sharp contrast with what one 
encounters in kaon 4-body decays [20l - [22| at least 
in two aspects: (i) in Ke4 decays the energy range 
is always far below resonance thresholds. As a re- 
sult, Bremsstrahlung is overwhelmingly dominant; 
and (ii) additionally, the kaon direct emission contri- 
bution can be estimated in chiral perturbation the- 
ory. 

In the present case, the phase space is larger 
and decays like D° — > V°£ + £~ have nonnegligi- 
ble branching ratios [23J. As a result, D° —> 
V®(h\h 2 )i + i~ are expected to play a significant 
role. However, unlike kaon decays, such contribu- 
tions cannot be estimated with effective field theories 
and one has to resort to hadronic models. Therefore, 
one has to exercise some caution when interpreting 
the results: hadronic effects are genuinely nonper- 
turbative and as such uncertainties are difficult to 
estimate. One should bear in mind that the results 
of this Section are not on the same solid ground as 
the Bremsstrahlung contributions discussed above. 

Quite generally, we will be considering underlying 
processes of the form depicted in Fig. [2] where 

T^(p,q)=i 3 J d 4 xe^ x+lc i-y\o\r c (o)jV( x )jj(y)\o) 

4-V 1 iV I 4.V 1 J.V 

= h q^lu + t 2 q^Pu + t 3 + t 4 e^xpPxQp 
+ tJp l iPv + tXp l ,qu (10) 




FIG. 2: Diagrammatic representation of the hadronic 
model adopted, with the weak T M „ (p, q) and strong 
B l _ l (p\,p'2) correlators. 



and 

B a (p U p 2 ) = i 3 Jd 4 xe^ x+ ^-y\0\jr(0)Ji(x)J 2 (y)\0) 

= b\pi a ~ b 2 p 2a (11) 

parameterize the most general tensorial decomposi- 
tion of the D° -> V{p)j(q) andV(p) -> h 1 (p 1 )h 2 (p2) 
vertices, tf and are functions of the kinemati- 
cal invariants of the problem, namely (q 2 ,p 2 ) and 
bj(p 2 ), where the upper index denotes their depen- 
dence on the exchanged resonance. 

The hadronic vector can be expressed as 

H»( Pl ,p 2 ,q) = ( hl h 2 \H\V)^%(V<y*\H\D°) 

Pv\V ) 

(12) 

where we will describe particle widths with a Breit- 
Wigner propagator, i.e., Pj(q 2 ) = q 2 — m 2 + iTjirij, 
and 

(h 1 h 2 \n\V)=B^ Pl ,p 2 )e^Hp) 
(Vr\n\D°) = T^(p,q)e<P'(p)ep>{q) (13) 
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In order to evaluate the matrix elements above we 
need to determine the functions bi and U. The for- 
mer can be easily determined from the experimental 
V — s- h\hi decays. Notice that for the equal mass 
case, charge conjugation invariance (in our case V is 
neutral) imposes that b\ = b\ . When ^ nih. 2 
the form factors differ by a term proportional to the 
mass difference (see Appendix [Bj . In general, it is a 
good approximation to neglect this effect and hence- 
forth we will use that b\ = b\ = by for all decay 
modes. We will also assume that the momentum de- 
pendence is soft and that to a good approximation 
b v (p 2 ) = 6y( m v0 = b v . 

The V — > h\h,2 decay width is 

r = J-b 2 v m 5 v \ 3 / 2 (m 2 v ,m 2 hl ,m 2 h2 ) (14) 

457T 

and comparison with the experimental determina- 
tions [6f gives 



b p = 5.92 GeV; b K . 



5.46 GeV; b d 



4.41 GeV 
(15) 



The contribution of the w(782) to D° -> n + TT^£ + £" 
is extremely suppressed (bu/b p ~ 3%) and will be 
neglected. 

The evaluation of the weak vertex requires the 
Ac = 1 effective hamiltonian. For the decay chan- 
nels that we are considering, the relevant opera- 
tors for the Cabibbo-allowed (C), singly Cabibbo- 
suppressed (SCS) and doubly Cabibbo-suppressed 
(DCS) transitions are 



ltA,-l — 7= 



rt/\r—l — — 7= 





Gf 


Ac=l — 


7f ■ 


SCS 




Ac=l — 


V2 i 


DCS 


G F r 


Ac=l — 


V2 L 



XdsG 



(da) 



ids 



where Xj = V c *V uj , X sd = V* s V udl X ds 



) sd = {vrf ll c) L (sYd)L 
Qd = (u-y,,c) L (d^d)L 
Q s = (u r y fi c) L (sj' J 's) L 



(16) 
V* d V U8 and 






1 0.1 



(17) 



In the previous equations we have neglected the pen- 
guin operators. To proceed further one needs to 
make some approximations. In this work we will 
assume that the factorizable contributions are the 
dominant ones, which means that 



(vr\j? 3 4 c \D n ) = (v\Jt 3 \o)(r\4 c \D o 
+ ( 1 *\jm)(v\jr^ 



(18) 



FIG. 3: Long-distance contributions to the differen- 
tial decay width (in arbitrary units) in the (mf;, m|^) 
plane for the different decay modes (from top to bot- 
tom: K^ir , -k + ty~ and K + K~ , respectively), ran and 
m,hh ire given in GeV. Above we show the e + e~ modes. 
For the dimuon case the only difference happens for the 
Bremsstrahlung , which is strongly suppressed due to lep- 
ton threshold effects. 



6 



where ij 
that 



dd,ss,sd,ds. We will further assume the ti(q 2 ,p 2 ) form factors is 



(i) The weak annihilation contribution (third 
line) is negligible compared to the spectator 
ones (first two lines) for all the processes to be 
considered (see, e.g., [24{ for a detailed discus- 
sion). 

(ii) The Zweig rule is at work, i.e., flavor annihila- 
tion is suppressed and a possible enhancement 
due to final state interactions is excluded. 

(iii) The photon is created mainly through vec- 
tor meson exchange and we will neglect a di- 
rect photon coupling. We will consider the 
exchange of the lowest-lying neutral states 
(p,u),<f>) to be dominant. 

The last point implies that 



Pv>{q 2 ) 



x (V'V\JUT\D°) 



(19) 



and requires the electromagnetic couplings of vector 
mesons, which can be inferred from 

n Vl = -\(v,„vn + J£- F ^(Qvn (20) 

where fv is defined by 

(V(k,e)\J^\0) =f v m v e;(k) (21) 

For phenomenological purposes we will break the 
SU(3) symmetry and take the experimental values 
for the decay couplings. The interaction term in the 
previous Hamiltonian will thus be replaced by 



%™ = -e 



rrir 



3m u 



V2U ( 
3m A 



(22) 



To complete the picture, we need the 
(V(p)\J^ c \D°(P)) matrix elements, which can 
be parameterized as (p + = P + p, k = P — p) 

(V(p,e)\jf\D\P)) = £> 1 (fe 2 ) P+M + D 2 (k 2 )k^ 

+ D 3 (k 2 )e* p + iD 4 (k 2 )e p „ XpP » + k x e<>* (23) 

A determination of the different form factors Di(k 2 ) 
is given in Ref . [25| . We refer to Appendix [B] for 
details. Matching Eqs. ([H]) and dT3j the result for 



m D + m p 
J y (q 2 )i 2 (p 2 )+^M/(<Z 2 )i 2 (g 2 ) 



t3(q 2 ,p 2 ) = i£,v(m D + nip) 



tl{q\p 2 ) 



J y (< Z 2 )i 1 (p 2 )+^T^( (Z 2 )i 1 (g 2 ) 
26/ 



m D + m p 

jVf n 2 



J v (q")V(p 2 )+r, v W(q')V(q 2 ) (24) 



where = Oa^^f") V V = 0,1 depending on the 
channel and 



J V (q 2 ) 
W(q 2 ) 
A 2 (k 2 ) 
il(fc 2 ) 
V(k 2 ) 



= q 



m p- p p(<7 2 ) 3m^P ti; (g 2 ) 



fvmv 



= q 



' ft , f 2 m ' 

P p {q 2 ) iP^{q 2 ) 3P (g 2 ) 



hA2m 2 



,12 



L A2 



h A im 2 A1 
-k 2 



'Al 



h vi m 2 /1 
k 2 



(25) 



l Vl 



As side remarks, we note that: (i) the first and sec- 
ond terms in Eq. (IM1) correspond to the first and 
second lines of Eq. (fT5)) . respectively. The absence 
of ^-exchange in J v (q 2 ) is thus a consequence of the 
Zweig rule; and (ii) the global prefactor (mo + m p ) 
in Eqs. (|2"4"|) results from considering the p(770) and 
cj(782) nearly degenerate. 

Bringing all the pieces together and matching to 
Eq. (H|) the form factors take the form: 



F, 



(V) 



.02\_q-Pi+aX 2 q-P2 + aX l 



Pv( P 2 ) 



rpty) _ .aJiQ ■ Pi + a 22q ■ P2 + 0-32 
r 2 —i 



PviP 2 ) 



(v) = q 4 

3 Pv(p 2 ) 



(26) 



where 



O&feV) = ~4 2 (q 2 ,p 2 ) = -ib V t%(q 2 ,p 2 ) 
a^(q 2 ,p 2 ) = ~al 2 {q 2 ,p 2 ) = - t b V t^(q 2 ,p 2 ) 
aUq 2 y)=-2b v t][(q 2 y) (27) 
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Results and discussion 



tors get therefore reduced to 



In order to proceed to a numerical analysis for the 
different decay channels we will take 26] 



2i£ 



K * 



fp = 216 MeV, 
/„ = 187 MeV, 



f K , = 220 MeV, 

U = 215 MeV (28) 



for the vector meson coupling constants together 
with Eq. ([15)1 for the V — > h^h^ couplings. In or- 
der to determine the parameters entering the form 
factors in Eq. (1251) , experimentally one has access to 
the residues [a, l23] : r v = hyi/h/a = 1.7 ± 0.8 and 
f'2 = hA2/h>Ai = 0.9 ± 0.4. Based on this results we 
will take hvi = 0.94 and hju = 0.55 and consider the 
couplings to be flavor-blind. Regarding the poles, 
we will adopt the usual values, namely mvi = 2110 
MeV and m A1 = m A2 = 2530 MeV [Hj]. We wil1 
also assume that the Wilson coefficient C2 — —0.55 
for all the channels |24j |. The remaining parame- 
ters (CKM matrix elements, decay widths, etc.) are 
taken from Q. 

The main results are collected in Table Q] and 
Fig. [3J In Table Q] we have listed the contributions 
of Bremsstrahlung, electric and magnetic emission 
to the branching ratios for the different channels 
(the electric-Bremsstrahlung interference is only at 
the 1 — 2% level and has been omitted). The am- 
plitudes turn out to be largely dominated by the 
electric pieces (Bremsstrahlung and electric emis- 
sion for the e + e~ decays, or only electric emission 
for the decays) with weights that depend on 

the channel under study. Fig. [3J shows the dif- 
ferential decay widths for the different channels in 
the (to^to, 2 ^) pl ane - Analytical expressions thereof 
can be found in Appendix [Cj Notice that, in agree- 
ment with Low's theorem, the low-q 2 region is dom- 
inated by the Bremsstrahlung. The resonance con- 
tributions (or any other contributions) should have 
a smooth low-g 2 behavior. Therefore, we do not 
quite understand the low-q 2 enhancement found in 
Ref. [H for D° -> Vt+tr, but it definitely cannot 
be ascribed to factorization, as suggested in (29j . 

In the following we will discuss each separate chan- 
nel in turn, focussing on their specific features. 



D 



K-n+l+r 



The Bremsstrahlung contribution is listed in Ap- 
pendix [C] The dominant resonance contribution 
to this Cabibbo-allowed decay comes from D° —> 
K*V', with subsequent decays K* — > K~tt + and 
V' — s- 7* — > . One can easily get convinced that 
strangeness conservation forbids the second line in 
the factorization formula of Eq. (|18 j) . The form fac- 



m D + m p 

tf (q 2 ,p 2 ) = ifr. (m D + m p )J K ' (q 2 )!^ 2 ) 
2£,K* tK* fjt\-Crr„ 2 



m D + m p 



(29) 



The combination of long-distance contributions is 
shown in the top panel of Fig. [31 where one can see 
the hadronic emission sitting (mostly) at (to 2 , m 2 K »). 
The w(782) contributions is also present but too 
small to be noticed by the naked eye. As expected, 
the Bremsstrahlung and resonance regions lie rather 
far apart, which explains why their interference is 
negligible. For the e + e~ case, the total branch- 
ing ratio adds up to 1.6 • 10 -5 , with a slight domi- 
nance of the Bremsstrahlung, while for we find 
6.7 • 10~ 6 , coming largely from the electric emission. 

The same qualitative features apply for the doubly 
Cabibbo-suppressed modes. In this case, branching 
ratios add to 5.5 • 10~ 8 for electrons and 1.8 • 10~ 8 
for muons. 



2. D° -»■ TT + TV-l + r 

The direct emission contribution for the D° — > 
7r + 7r~7* decay comes mainly from the near-resonant 
decays D° —> p°j* and D° — > ujj* . As we argued in 
the previous subsection, for phenomenological pur- 
poses it is a good approximation to neglect the u> 
contribution due to its tiny branching ratio into two 
pions. 

Fig. 131 shows both Bremsstrahlung and direct 
emission contributions. The resonance region now 
shows distinct peaks at (to 2 , to 2 ) and (to 2 ;, to 2 ) of 
comparable size. Contrary to the Cabibbo-allowed 
mode, in this case both spectator pieces contribute 
and Eq. ([Mf assumes its more general form for the 
t p 

The branching ratios for the e + e~ and are 
respectively 2.0 • 10~ 6 and 1.5 • 10~ 6 . We note that 
our results for the direct emission are in very good 
agreement with Ref. [2{|, where the 3-body decays 
D — > pl + l~ were considered. Our results are espe- 
cially interesting in view of the prospects of LHCb 
to be able to detect signals at 10 -6 in Cabibbo- 
suppressed D° — > h±h,2 /i + decays [ill ]. 



3. D° -> K+K-l+r 

The direct emission contribution comes in this 
case mainly from a near-resonant D° — > 07*. Simi- 
larly to the Cabibbo-allowed case, only the first line 
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survives in factorization and therefore the form fac- 
tors reduce to Eq. (|29l) . with the obvious substitu- 
tion K* — !• (j). 

The resonance region is peaked at {m 2 p , m^), such 
that Bremsstrahlung and direct emission populate 
opposite sides of the phase space in the Dalitz plot. 
The total branching ratios for the e + e~ and fi + fi" 
are respectively 6.5- 10 -7 and 1.1 • 10~ 7 , with a slight 
dominance of the Bremsstrahlung for the electron 
case. 



III. SHORT-DISTANCE EFFECTIVE 
HAMILTONIAN 

In the previous Sections we have studied the (dom- 
inant) long-distance contributions to the D° —> 
hfh,2£ + £~ decays. In this Section we turn our atten- 
tion to short distances. Our main aim is to single out 
an observable without long-distance background. In 
charm decays, such an object is automatically a new 
physics probe, since the standard model contribu- 
tion is not only loop-suppressed, but additionally 
shrinked by the GIM mechanism. 

The c — > ul + l~ transition inside the standard 
model is mediated by electromagnetic and Z pen- 
guin diagrams and by W box diagrams at one-loop 
order. At energies right below the charm threshold 
it is described by the following effective lagrangian 

^.(c^ui+r) _/w(c-wy) ^ F \^ \ \n^n x^'n 

tt eff — H ef f -j= } j Aj |Q 9 ^9 + ^10^10 

v i 

(30) 

where = V^Vui, i — d,s run through the down- 
type quarks and 

l07T^ 

Oio = ^iu^P L c)(I^ 5 £) (31) 

The full list of operators in Ji^^" 1 ^ can be found 
in 23]. The leading contribution in c — > try comes 
from the electromagnetic penguin operator Q-jl 

Qil = i^mX(ua^P R c){l^e) (32) 

The remaining operators are 0(a s )-suppressed [30l ] 
and will be neglected. We have also omitted Qtr 
and Q' g 10 , which are weighted by m u and can be 
safely neglected. 

The RG running of the previous operators was 
first computed in [31( to the leading logarithm ap- 
proximation. QCD corrections to c — > 1*7 were 
shown to lead to a significant enhancement of 



CVl [H[ , where the strongest effect comes from the 
two-loop contribution: mixing with mainly Q 2 soft- 
ens the GIM suppression from power-like to loga- 
rithmic [32|- This phenomenon was first pointed 
out in [33| and later applied to B decays [34| and 
kaon decays [H, [36j]. Despite this enhancement, 
when it comes to c — > ul + l~, Qjl is overshad- 
owed by Qg at least by one order of magnitude, 
since its mixing with Q 2 happens already at tree 
level ■ Although one-loop QCD corrections are 
significant (3?J and effectively reduce the value of 
C9, it still constitutes the dominant contribution 
to the short-distance standard model estimate of 
c — > ul + l~ 38]. In contrast, Q10 does not mix with 
Q2 and is invariant under the renormalization flow, 
which makes its value extremely suppressed in the 
standard model: 

C 10 (m c ) = C w (m w ) - —§- (33) 
m w 

For inclusive c — > ul + l~ transitions, the branching 
ratio has recently been estimated at 3.7 • 10~ 9 for 
electron-positron dilepton pair [38| . For the dimuon 
case, it is resonable to expect a factor 5 suppres- 
sion. Taking into account our results in table HI this 
entails that short-distance effects have a negligible 
impact on the decay width, as expected. Even if 
one is extremely conservative and makes the com- 
parison only with the Bremsstrahlung, short dis- 
tances might be competitive only for the dimuon 
Cabibbo-suppressed modes. However, disentangling 
short from long distances in those cases is extremely 
challenging: since the bulk of the short distances 
comes from Cg, the short-distance piece has the same 
distribution in the Dalitz plot as the long-distance 
one, i.e., the short-distance contribution will simply 
pile up on top of the long-distance contribution. 

IV. ANGULAR ASYMMETRIES 

The most general angular distribution for a 4- 
body decay can be parametrized, using the Cabibbo- 
Maksymowicz set of variables, as in Eq. (jTJ]), which 
we repeat here for convenience: 

d 5 T 

= Ai(x) + A 2 {x)sj + A^{x)s t c (tl + Ai{x)s 2 pc t)> 

+ A 5 (x)s e c cf> + Ae(x)ce + A 7 (x)siS cf> 

+ A 8 (x)s 2 £s c f > + Ag(x)sjs 2< t, (34) 

The pieces contributing to the decay width were 
studied in Section UU and correspond to the first 
line above. The remaining angular structures can 
be generated by interference effects, either at the 
hadronic or at the leptonic vertex, and can be probed 
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Mode 




(e+e- 









A* 






3.9 


io- 7 (^ 


-2%) 


3.4 


10 




2%) 


7T + 7T~ 


1.3 




-6%) 


1.1 


10 




8%) 


K+K~ 


7.9 


10~ 9 (. 


- 1%) 


6.9 


10 




6%) 


K+TT~ 


1.1 


10" 9 (r 


-2%) 


9.6 ■ 


10" 




5%) 



TABLE II: Branching ratios for the interference between the electric and magnetic terms assuming 8w — 7r/4 for the 
weak phase. In parenthesis we show the value for the T-odd asymmetry A$. 



with different angular asymmetries. It is important 
to emphasize that, due to the rich angular struc- 
ture of Eq. (|34l) . angular asymmetries are not re- 
stricted to charge asymmetries. Moreover, there 
are asymmetries that, due to tiny standard model 
backgrounds, turn out to be clean probes of new 
physics. In this Section we will concentrate on two 
such asymmetries. Far from being exhaustive, we 
just want to provide some exploratory indications of 
what are their expected signals in the Dalitz plot and 
parametrize its size in terms of a few short-distance 
parameters. 



A. T-odd asymmetry 



It was long noted, in the context of Kl — » 
7r + 7r~e + e~, that the interference between the mag- 
netic and Bremsstrahlung contributions can lead to 
large CP violation [39] . This interference is genuine 
of 4-body decays and contributes to As.g m Eq. (j3~4")l . 
Therefore, it can be singled out by an angular asym- 
metry in the diplane angle cf>: 



1 f 27T dT 

A* = <sgn( W )> = - J - 



d<t>* (35) 



where we have defined the piece-wise angular inte- 
gration 
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tt/2 



tt/2 



3tt/2 



2tt 



>3n/2 

(36) 

This asymmetry actually collects not only the 
Bremsstrahlung vs. magnetic interference, but also 
the electric vs. magnetic one. In kaon physics the 
former dominates and, since the Bremsstrahlung and 
magnetic contributions have different strong phases, 
A^ becomes a probe of long-distance physics (see, 
e.g., H). 

In D decays, in contrast, the electric vs magnetic 
interference is expected to be dominant, because 
both pieces have the same structure in the Dalitz 
plot. In comparison, the Bremsstrahlung vs mag- 
netic interference here is severely suppressed, even 
allowing for large strong phases. Interestingly, the 



electric vs magnetic interference can only be nonzero 
in the presence of weak phases. In Section|TT]we eval- 
uated the hadronic matrix elements for the dominant 
Qi four-quark operator. Weak phases emerge when 
other four-quark operators are also considered, such 
as Qi or the QCD penguins, so that their presence is 
naturally expected. In the standard model, however, 
they are extremely suppressed and can only become 
sizeable in the presence of new physics. Therefore, as 
opposed to kaon physics, in D physics A^ is a probe 
of new physics, with signals mostly concentrated on 
the resonance region of the Dalitz plot. 

In Fig. @] we show the differential electric-magnetic 



interference as a function of 



for the different 



decay modes. For comparison, we have included 
the long-distance background (dashed line) . For the 
sake of illustration we have picked Sw ~ 7r/4 as a ref- 
erence value for the weak phases. As expected, the 
contributions are sizeable close to the exchanged- 
resonance peaks, which vary depending on the fi- 
nal dihadron state. The interference between mag- 
netic and Bremsstrahlung (not shown in the plot) 
is only appreciable at the dilepton threshold and 
it amounts to only 1% of the integrated asymme- 
try. In Table |TT] we have listed the resulting val- 
ues for A^ for 5w ~ tt/4. Notice that in that case 
the asymmetry hovers in the 0(1 — 10%) window, 
depending on the channel, with D° — > ir + ir~ be- 
ing the most favored decay mode. Quite generically, 
dimuon modes give bigger signals than e + e" modes, 
which should be taken as an additional motivation 
to study the dimuon decays by LHCb, specifically 
D° — > 7T+ ir~ fi + fi~ . 



B. Forward-Backward asymmetry 

In the previous subsection P violation was induced 
in the hadronic vertex. It is also interesting to con- 
sider P violation in the leptonic vertex. 

Consider the matrix element stemming from the 
scmilcptonic penguins Qg and Qw, namely 

Mf D =^{k + ,k^)H ll {p l ,P2,q) 
MfJ =Z w L» 5 (k + ,k_)H^ Pl ,p 2 ,q) (37) 
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FIG. 4: The angular asymmetry (see Eq. A35\) ). The solid lines show the normalized differential distribution as a 
function of (in GeV 2 ) for the different decay modes. The dashed line corresponds to the differential decay width, 
which is included for comparison. From top to bottom: K~-k + , tt + -k~ , K + K~ and K + tt~ modes. Left and right 
panels collect, respectively, the e + e~ and /i + fi~ modes. 
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where 

L M (fc+,fc_) = u{k^v{k+) 
L^ 5 (fc+,fc_) = u(/fc_)7^7 5 i>(fc + ) (38) 
The short-distance hadronic tensor is defined as 

n^ PuP2 ,q) = (h 1 h 2 \j; c \D a ) (39) 

and 



£(9,io) — ~/=^6C(9,10) 



16tt 2 



Gf 

71 



AbC (9 , 10) (40) 



Above, the CKM unitarity relation A^ + A s + A& = 
has been used. The short-distance hadronic vector 
T~L[i admits a decomposition in terms of form factors 
•Fl.2,3 akin to the long-distance one in Eq. (@}. As 
in the long-distance analysis of Section UH the dom- 
inant contribution comes from the near-resonant re- 
gion. Thus, 

(h^J^D ) = ( hlh2 \n\V)p^{V\4 c \D°) 

(41) 



Following the steps of Section ITTl and the expressions 
given in Appendix [B] the form factors can be shown 
to be: 



v 92(q 2 )q ■ (pi -P2) +g3(g 2 ) 
Pv(p 2 ) 

,v92(i 2 )q ■ (pi -P2) - fa{q 2 ) 



Pv{p 2 



-1b x 



Pv(p 2 



(42) 



where 



92 (q 2 



2i 



-A 2 (q 2 ) 



93(l 2 ) = i{ m D + m p )Ai (q 2 ) 
g4q 2 ) = V{q 2 ) 



(43) 



m D + m f 

and we have used that b\ = b\ = b v . 

The interference with the long-distance (photon- 
mcdiatcd) contribution studied in Section HTl reads 

Re[Ml D M<£$] = ^ £ Gi&mFiFj-FjFt] (44) 

where Gij are given by [22j 
G12 = e uv \ p PiP2Q x q p 



12 
G13 



- [pi ■ Q(q 2 Pi -P2~p\- qp2 ■ q) 

+P2 ■ Q((P1 ■ qf - m 2 hl q 2 )] 

G 23 = [p2 ■ Q(q 2 P2 -p\-P2- qpi ■ q) 

+Pl ■ Q((P2 ■ qf - m 2 h2 q 2 )] 



(45) 



The term proportional to G\ 2 describes the inter- 
ference between the electric components of long 
and short distances, while G13 and G23 collect the 
electric vs magnetic interference terms. They con- 
tribute to the second line in Eq. (|34l) . The forward- 
backward asymmetry in Of. 



A 



FB 



(sgn(Q)) 



d y 



, dr 

dy- 

-1 



dy 



(46) 



singles out Aq and it is sensitive to the electric vs 
magnetic interference. Due to the tiny value for C10 
in the standard model discussed above, a nonvanish- 
ing Aps is a rather clean test for new physics. 

In Fig. [S] we show the differential electric-magnetic 
interference as a function mf t for the different de- 
cay modes. The figure actually shows the absolute 
value of the interference: since strong phases are 
accounted for as Breit-Wigner widths, the contribu- 
tion flips sign at each resonance pole. For compari- 
son, we have included the long-distance background 
(dashed lines). Notice that the angular asymme- 
try stays rather constant from the dilepton thresh- 
old up until the resonance peak, beyond which it 
falls off rather steeply. In order to be quantitative 
we have considered a new physics scenario in which 
C10 ~ 0(1). This is the case, for instance, in the 
Little Higgs model studied in Refs. [HI, The 
expected size of Afb for the different decay modes 
is summarized in the first two columns of Table IIIII 
Notice that since C10 is flavor-universal, i.e. not 
subject to the CKM hierarchy, the signal steadily 
increases from the Cabibbo-allowed to the doubly 
Cabibbo-suppressed modes, where it can reach the 
10% level for the dimuon mode. 

This last observation suggests to consider also the 
interference between Cg and C10 in new physics sce- 
narios where Cg ~ C10 ~ 0(1). In this case the 
CKM hierarchy is absent altogether: both Cg and 
C10 are flavor-blind. The results are summarized in 
Fig. [5] and the second column of Tabic IIIII Analyti- 
cal expressions can be found in Appendix [C] Notice, 
as compared to the previous case, that the gradual 
dominance as one goes from the Cabibbo-allowed to 
the doubly Cabibbo-suppressed is more pronounced, 
as expected. However, we note that even with this 
huge enhancement in Cg and C10, the purely short- 
distance interference is only appreciable for the dou- 
bly Cabibbo-supressed modes. 



CONCLUSIONS 



The recent evidence of CP violation in D° — > 
h + h~ decays well above the expected standard 
model prediction makes the study of D decays a pri- 
ority in the search of new physics. In this paper we 
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Mode (e + e )sd-ld )sd-ld 


(e + e )sd-sd )sd~sd 


K--JT+ 4.5 ■ 10" 8 (~ 0.3%) 4.0 • 10" 8 (~ 0.2%) 
ty+%- 2.8 ■ 10~ 8 (~ 1.4%) 2.6 • 10" 8 (~ 1.8%) 
K+K- 2.8 ■ 10" 9 (~ 0.4%) 2.4 • 10~ 9 (~ 2.1%) 
K + ir~ 2.4 • 10 -9 (~ 4.3%) 2.1 • 10" 9 (~ 11.4%) 


1.7 • 10" 9 (~ 0.01%) 1.4 ■ 10" 9 (~ 0.009%) 
3.0 • HT 9 (~ 0.1%) 2.6 • 10~ 9 (~ 0.2%) 

4.4 ■ 10- 10 (~ 0.07%) 3.6 ■ 10- 10 (~ 0.3%) 
1.7 • 10" 9 (~ 3.0%) 1.4 • 10" 9 (~ 7.8%) 



TABLE III: Branching ratios for: (i) the interference between short distances (Qw) and long-distance direct emission 
(denoted SD — LD) and (ii) the pure short- distance interference between Qg and Qio (denoted SD — SD) for the 
different decay modes. In parenthesis we include the value of the forward-backward asymmetry Afb- As reference 
values for the penguin coefficients we have chosen Cg t io ~ 4 (see Ref. fluj). 



have provided the first detailed analysis of the rare 4- 
body decays D° -> h+h^i+i^, (h = n,K;£ — e,(i) 
in the standard model. We have studied the dom- 
inant long-distance contributions (Bremsstrahlung 
and hadronic effects) in the (m^m^) Dalitz plots 
and the total branching ratios, which turn out to be 
substantially larger than previously estimated. Both 
the Dalitz plots and the associated branching ratios 
should be useful tools in view of the upcoming anal- 
yses of LHCb and BESIII. 

For the Cabibbo-allowed, singly Cabibbo- 
suppressed and doubly Cabibbo-suppressed modes 
one finds 

Br[D° -> K--K + l + r] ~ 1CT 5 
Br[D° -> Tr+w-e+J?-] ~ 10~ 6 
Br[D° -> K + K~£ + £~] ~ 10~ 7 
Br[D° -> K + tt-£+£-} ~ 10" 8 (47) 

where an important contribution comes from the 
near-resonant processes D° —> V£Vf, with subse- 
quent decays VJ* hih 2 and V e * -t 7* ->• £ + £~. As- 
suming factorization, lowest meson dominance and 
strong phases coming mainly from the resonance 
widths, the different form factors involved can be 
determined from experimental input. Our results 
for the hadronic contribution are in agreement with 
Ref. [29j | . which is quite a nontrivial agreement, given 
that the assumptions and methods going into both 
analyses are rather different. Finally, we want to em- 
phasize that our results comply with Low's theorem, 
i.e., the main contribution in the low dilepton invari- 
ant mass region comes from the Bremsstrahlung. 

After having determined the standard-model con- 
tribution to D° —> h\lru l £ + £~ we have also ex- 
plored signals for new physics detection. In par- 
ticular, we have shown that two angular asymme- 
tries, namely the T-odd diplane asymmetry and the 
forward-backward dilepton asymmetry can provide 
direct tests of new physics due to tiny standard 
model backgrounds. Using some sample values for 
the (unknown) short-distance parameters, we show 
that effects can generically reach the 10% level. 
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Appendix A: Kinematics 



We define p = p\ + p 2 and q = k + + fc_ as the 
momenta of the dihadron and dilepton pairs, respec- 
tively. Then one can obtain for the phase space 




Above we have defined = A(p 2 , m^, m^ 2 ) and 
Ad = \(p 2 ,m 2 D ,q 2 ) where A(a, 6, c) = a 2 + b 2 + c 2 - 
2ab — 2ac — 2bc. The angular variables are denned 
as in Ref. [17]: if Pi is the hi momentum in the 
dihadron CM system; k + the £ + momentum in the 
dilepton CM system; n the direction of the dihadron 
system as seen from the D° rest frame; and p± and 
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FIG. 5: The forward-backward asymmetry Afb (see Eq. ^46\ )). The solid lines show the normalized differential 
distribution as a function of va\ (in GeV 2 ) for the different decay modes. The spiky line corresponds to the absolute 
value of the long vs short-distance interference (with Cio = 4j, while the smooth solid line is the pure short-distance 
interference (with Co, = Cio = A), The dashed line corresponds to the differential decay width, which is included for 
comparison. From top to bottom: K~ir + , tx + -k~ , K + K~ and K + tt~ modes. Left and right panels collect, respectively, 
the e + e~ and fi + fj,~ modes. 
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k:J: the components of pi and k + perpendicular to 
n, then 



a n P 1 a 

cos Oh = — — ; cosw 

IPil 



fi • k_i 



|k+| ' r \vi\H\ 
(A3) 

Intuitively, 9e is the angle between the £ + momen- 
tum and the dihadron system as measured from the 
dilepton CM while <f> is the angle between the di- 
hadron and dilepton planes. 

The final result for the phase space therefore reads 



rf 5 $ = 



2 14 7r 6 m 2 p 2 



4mj 



x dp 2 dq 2 d cos 9 hd cos 9 edcj) (A4) 
where the range of the kinematical variables is 
4m| < q 2 < (m D - {mux + m h2 )) 2 
(m h x + m h2 ) 2 < p 2 < (m D - ^/q 2 ) 2 
< (9 h , 6 t ) < it 

< </> < 2tt (A5) 

The relevant kinematic products can be expressed 
as: 

1/2 2 2 \ 

Pl-P2 = -(p -m hl -m h2 ) 

Q-Pi,2 = \{m 2 D -p 2 - q 2 )C± T -TzV^h^D cos 9 h 
4 4p z 



Q ■ Pl,2 = Pi 



An— cos9p ± 



2 2 2 

m D — p — q~ 



p A 



■ cos Vh cos 9 1 =F 



2p 



j- \J q 2 p 2 sin 9h sin 9e cos < 



» , i nP _ p 2 - ( m li - m D 



^vXpPxP^q Q 



4p 4 

x fti sin 9h sin (9^ sin < 



D q 2 p 2 



where 



ft 



and x 2 = — m h2- 




(A6) 
(A7) 



The most general parametrization for the remaining 
matrix elements is 

(V(p, e) | Jf | D° (P)) = Dx (k 2 ) P+p + D 2 (k 2 ) k, 

+ D 3 (k 2 )e; + iD±(k 2 )e^ p p\k x eP* (B2) 

with p + = P + p and k = P — p. A similar ex- 
pression holds for (V'\J™ C \D°). The divergence of 
the current is proportional to the difference of quark 
masses, which can be parametrized as the squared 
mass difference of the hadrons. In other words, 



k»(V(p,e)\j; c \D Q (P)) ~ (m% - m 2 v )k ■ e* 
which is satisfied if 

Dx{k 2 ) = Ai(fc 2 )fc-e* 



(B3) 



D 2 {k 2 ) = A 2 (k 2 )k ■ e* 



D 3 (k 2 ) = A 3 (k 2 )(m 



D 



k 2 

m 2 v ) 



(B4) 



In particular, the second equation above implies that 
^2(0) = 0, such that the divergence is avoided. If 
one now shifts A 2 — > A 2 — A 3 — A\ one obtains 

(V(p)\4 c \D°(P))=A 1 (k 2 )k-e* 



a _ (m 2 D - ml) 
+ k 2 



+ A 3 (k 2 )(m 2 D - m 2 v ) 
+ iD 4 (k 2 )e^ xp p + „k X E 



k 2 



k-e* 



k 2 



(B5) 



which can be straightforwardly compared with the 
parametrization of Ref. 25]. The form factors are 
related as 

Ai(k 2 ) = -i(m D + my)" 1 A 2 {k 2 ) 

A 3 {k 2 ) = i(m D - m v )- 1 A l {k 2 ) 

D A {k 2 ) = i(m D +m v y 1 V(k 2 ) (B6) 

The hatted form factors are determined in Ref. [25[ 
in the nearest pole approximation (vector meson res- 
onance) as 



Appendix B: Hadronic parameterization 

The hadronic matrix element for the weak vertex 
can be written, using vector meson dominance and 
assuming factorization, as 

(V(pmq)\J$J?\D a (P)) = Y.mv^v)^^ 

V V'\H ) 

x {fvmvefciV'lJFW + fv'mv'^WJFlD )} 

(Bl) 



V(k 2 ) 



h vi m 2 vi 
k 2 



'VI 



A 2 (k 2 



Ax(k 2 ) : 

hA 2 m 2 



hAxm 2 A1 
k 2 ' 



"AX 



A2 



(B7) 



L A2 



On the other hand, one has that 



(V(ph*(q)\J^ c \D^P)) = T^(p,h)el{p)e*M 

= t v 2 (q ■ e* v )(p ■ e ;) + tj(e* v ■ C,) + t v 4 e^tl*el*p x q p 

(B8) 
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Comparison between Eqs. (IB1|) and (|B8I) yields the 
expressions appearing in Eq. (f24f in the main text. 

Regarding the strong vertex, the most general 
parametrization is 



B»{p + ,p.) = Bi(p 2 + K + B 2 (P 2 + V- 
where p± = p\±p 2 . Since 

P+»B> 1 {p + ,p-) ~ (m 2 hl - m 2 h2 ) 
this entails that 



(B9) 



(BIO) 



B l (pl) = K l {p\ 
B 2 (p 2 + ) = n 2 {p\ 



m hi - m 



1:2 



2 

P+ 



(Bll) 



where Ki(0) = to smoothen the divergence out. 
Therefore, 



B"(P1, P 2) = (K 2 + K lX 2 K 
EE b x p^ + b 2 P2 



(K2 



2 

(B12) 



where \ 2 



'hi 



m 2 l2 )P+ 2 ■ I n the equal dihadron 



mass case, b\ = b 2 , in agreement with charge conju- 
gation invariance. In the general case one finds that 
b\ ^ b 2 . However, since bi ~ O(GeV), the mass 
correction amounts at most to an 0(1%) correction 
and can be safely neglected. 



Appendix C: Differential decay widths 



cay widths for the generic case D° — > h^h 2 £ + £^, 
i.e., for rrihi ^ vtiy i2 - The singly Cabibbo-suppressed 
case, i.e., m^i — m} l2 , can then be worked out as a 
particular case. We start with 



dT = 



1 



2rau 



\M\ 2 d$ 



(CI) 



where the differential phase space rf$ is given in Ap- 
pendix [2] and 



M =M LD +M i s ) D +M< s 1 $ 



(C2) 



whose definitions can be found in the main text. It 
is useful to define the following functions: 



h+-(q 2 , P 2 ) 
h-(q 2 ,p 2 ) 
h- + (q 2 , P 2 ) 
h ++ (q 2 ,p 2 ) 



(m 2 D 
(m 2 D 



+ q 2 -p 2 ) 
q 2 -p 2 ) 
q 2 +p 2 ) 



P 



(C3) 



and the mass combinations 



l hl 



L hl 



'I i.2 



'1,2 



(C4) 



In this appendix we will provide analytical ex- 
pressions for the angular-integrated differential de- 



For the long-distance contributions with symmetric 
angular integration one finds 



d 2 T 
dq 2 dp 2 



d 2 T 
dq 2 dp 2 



Cb 



x log 



7 2 - 4m 2 



hi 



q 2 - 4m 2 2 



2(rh 2 -p 2 ) + q 2 

(X 2 h— +p 2 h + _Y - \ h X D ' ( X -'//__ - [rh- I-' - XhXo + 2^JEp 2 h+J 

X A h 2 __-{p 2 h + _-^JBf 
X 4 h 2 __ - (p 2 h+- + VhJo) 2 



(C5) 



El. 



p 2 (p 2 - 2fh 2 )X 2 D + 4 X 4 (A D + 3qY)X D \\a v 21 \ 2 + 4h—(h' + _ - Am 2 D q 2 ){\ h + 3 X 4 ) 



x Re 



V V* 




a 21 a 31 


+ 4 { 



+ 4^ p 2 (p 2 - 2rh 2 )(\ D + UqY) + 4 X *(\ D + 3q 2 p 2 ) \a^\ 2 



(C6) 
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d 2 r 



dq 2 dp 2 
d 2 T 



Cm 



Mag. 



\Pv\ 2 



(C7) 



dq 2 dp 2 



Br. -El. 



8h-+Re 



+ WRe 



*31 



x 2h + . 



(-iy2 X 2 h_ + 





'<4i~ 




Re 


+ 4 




[Pv. 



\ p 2 sr io ■ - * 2 fr— - 

(C8) 



(-iy2 X 2 -/i ++ + 8m2. 









Re 


~ a 31~ 


)] 


[Pv. 





where 

Cb-- 



Cm 
Ceb 



In 2 , 2\o V^h^DP 2 , . . , 2 

— (2m, +g )ft j |7W D | 



a 



{\h\d? /2 



18432^3 ( 2 ™ 2 + g 2 )ft- 



v 3/2 



2m, +g )fo 



\Mi 



(C9) 



and we have defined .Md = 



We 



have included for completeness the Bremsstrahlung 
vs. electric interference, even though its effect is 
extremely suppressed. Expressions for the mode- 
dependent form factors a V i , a^i and a\ can be found 
in the main text. 



For the short-distance contribution entering 
one finds that 



d 2 T 



dq 2 dp 2 
d 2 T 



C 



em 



El. -Mag. 



IW)I : 



■ sin SwRe 



,v „v 



dq 2 dp 2 



= Cbm sin S s Re 



Br. -Mag. 



at 



Pv 



(CIO) 



x log 



p 2 h+- + x 2 h- 



p 2 h^ 



X 2 h — + \A/Ad 



^\ / %Ad + |m 2 ll (m|,/i__ -p 2 h+J) + q 2 {x A + m 2 D p 2 ) - m 2 h2 q 2 h- 

(Cll) 



m h i <-> m h2 ,x 2 <-» ~X 2 



r 



where 5\y and 6s are, respectively, weak and strong 
phases and 



a 



,EM 



Cbm 



23047r 5 mf, 



\' i l 2 \ 

ti\^(.e- P 2 ) 



q 2 p e 



v 3/2 



96tt 9 / 



P$^(e-P 2 )\M D \ (C12) 



2,™3 



m; 



D 



Again for completeness we have included the 
Bremsstrahlung vs magnetic interference term, de- 
spite being strongly suppressed. 

Finally, we list the expression for the contri- 
butions entering the forward-backward asymme- 
try Afb, namely the terms Re[M* LD Mg^] and 



Re[Mf*M^]: 
d 2 T 



dq 2 dp 2 



d 2 T 



dq 2 dp 2 
where 



SD-LD 



SD-SD 



AD 

^FB 

\Pv\ 2 



<r(2) 
\Pv\ 2 



Re 



Re 



V -> V* , '■V V* 
a 31< + a 31< 



"V ~v* 
a 31 a 4 



cW _ 



3/2 . 
P 4 



i 



ft \3/2 1 2 

(2) _ 4g4io R i x h a dQ 
-Pe 



(C13) 
(CM) 



(C15) 
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and we have used the short-hand notation Expressions for the short-distance coefficients £9,10 

and the short-distance form factors gi(q 2 ) can be 
"21 (<? 2 ) = - ibV (J2{q 2 ) found in the main text. 

« 3 W) = -ib V g 3 (q 2 ) 

al{q 2 ) = -2^5 4 (g 2 ) (C16) 
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